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ABSTRACT 

The  passivity  theorem  may  be  used  to  design  robust  controllers  for  structures 
with  positive  transfer  functions.  This  paper  extends  this  result  to  more  general  configu¬ 
rations  using  dissipative  system  theory.  A  stability  theorem  for  robust,  model-independent 
controllers  of  structures  which  lack  collocated  rate  sensors  and  actuators  is  given.  The  theory 
is  illustrated  for  non-square  systems  and  systems  with  displacement  senso^sy— : . . . 
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1  Introduction 


In  this  paper  we  are  concerned  with  control  of  systems  with  second-order  dynamics  modelled 
by  the  following  system  of  ordinary  differential  equations. 

Mx(t)  +  Dx(t)  T  Kx(t)  =  F-u(t)  (1) 

y(t)  =  cdx{t)  +  Cvx(t)  +  Cax(t )  (16) 

For  structures,  this  is  an  approximation  to  an  (infinite-dimensional)  partial  differential  equa¬ 
tion  model  and  this  system  of  equations  may  be  of  very  high  order,  although  the  structural 
matrices  M,  D  and  K  are  positive  definite  and  typically  sparse  and  symmetric. 

Control  theory  for  time-invariant  linear  systems  which  are  described  by  first-order 
dynamic  equations  has  been  well  established  for  decades,  and  many  control  software  tools 
are  available  today  for  systems  written  in  first-order  forms.  For  applications,  engineers  can 
simply  convert  whatever  models  they  have  to  the  first-order  forms  and  then  use  the  existing 
tools  to  design  the  controllers. 

For  second-order  systems,  transforming  to  first-order  form  not  only  increases  the 
dimension  of  the  problem,  but  also  destroys  the  sparsity  of  the  structural  matrices.  Not  only 
physical  insight,  but  computational  efficiency  is  often  lost  in  conversion  to  first-order  form. 
Existing  control  analysis  and  design  software  may  not  be  able  to  handle  such  a  large  system. 
For  example,  solving  a  1000- by- 1000  Riccati  equation  is  practically  impossible  with  today’s 
numerical  techniques. 

There  are  basically  two  ways  to  address  the  controller  design  problem  for  a  large 
scale  system.  One  way  is  to  minimize  the  dimension  of  the  system  model,  through  some 
model  reduction  technique.  The  reduced  model  is  used  in  the  controller  design  and  some 
robust  design  methodology  is  used  so  that  the  controller  stabilizes  the  original  high  order 
model.  This  is  the  approach  behind  most  Hqo  design  for  structures  [6]. 

Another  way  is  to  design  a  controller  which  is  independent  of  uncertainties  in 
the  system  model.  The  advantage  to  this  approach  is  that,  unlike  the  first  approach,  the 
accuracy  of  the  original  high-order  model  (1)  is  not  assumed,  nor  is  an  accurate  system 
identification  required  This  is  very  appealing,  since  at  the  current  time,  no  accurate  model 
of  structural  damping  exists,  and  the  existing  models  of  stiffness  are  only  accurate  at  low 
frequencies.  Furthermore,  damping  is  difficult  to  identify  experimentally. 

Structures  with  collocated  rate  sensors  and  force  actuators  ( C'v  =  F)  are  passive, 
t'.e.  have  positive  transfer  functions  [2].  The  passivity  theorem  [1,  3]  implies  that  any  positive 
controller  will  lead  to  a  stable  closed  loop  system.  This  fact  can  be  used  to  design  control 
systems  which  remain  stable  despite  large  modelling  errors.  Recently,  several  researchers 
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have  incorporated  this  approach  with  other  robustness  theorems  in  robust  controller  design 
for  systems  which  have  collocated  rate  sensors  and  force  actuators  but  whose  transfer  func¬ 
tions  may  be  non-positive  due  to  computation  delays,  actuator  dynamics  etc.  [7].  However, 
the  restriction  of  collocated  sensors  and  actuators  is  a  stringent  one,  and  excludes  many 
applications.  In  particular,  systems  with  displacement  sensors  and  systems  with  a  different 
number  of  inputs  and  outputs  cannot  be  handled  with  this  approach. 

Since  altering  the  configuration  of  sensors  and  actuators  does  not  alter  the  fact 
that  structures  with  positive  damping  dissipate  energy,  the  question  arises  as  to  whether 
the  passivity  theorem  can  be  generalized.  For  instance,  when  a  mass-spring-dashpot  is 
attached  to  any  mechanical  system,  including  flexible  space  structures,  the  damping  of  the 
system  is  almost  always  augmented  regardless  of  the  system  size.  The  parameters  of  the 
mass-spring-dashpot  are  arbitrary,  model-independent  and  thus  insensitive  to  the  system 
uncertainties.  Knowledge  of  the  system  model  may  be  used  in  adjusting  the  parameters  in 
order  to  satisfy  system  performance  requirements.  However,  changes  or  errors  in  the  system 
will  not  destabilize  the  system  because  it  is  an  energy-dissipative  device.  The  question 
arises  as  to  if  there  are  any  feedback  controller  designs  using  sensors  and  actuators  which 
behave  like  the  mass-spring-dashpot  device.  The  answer  is  positive  as  shown  in  [4],  wherein  a 
conventional  engineering  approach  was  used  to  design  active  controllers  which  behave  like  the 
mass-spring-dashpot  device.  Although  the  conventional  approach  provides  physical  insights, 
it  lacks  the  theoretical  basis  necessary  for  application  to  more  complicated  systems  such  as 
nonlinear  flexible  robots.  This  paper  is  motivated  by  the  desire  to  provide  a  rigorous  theory 
to  support  this  argument  based  on  physical  intuition. 

In  the  next  section  we  summarize  results  on  dissipative  system  theory  obtained 
mainly  by  Willems  [8,  9],  A  stability  theorem  which  generalizes  the  passivity  and  small 
gain  theorems  is  given.  In  order  to  facilitate  the  application  of  these  ideas  to  non-linear 
models  of  structural  dynamics,  no  assumption  of  linearity  is  made.  In  subsequent  sections 
we  apply  this  theory  to  derive  conditions  for  robust,  model  insensitive  control  of  (1)  non¬ 
square  systems  with  rate  sensors  (2)  systems  with  displacement  sensors.  As  an  example,  we 
present  a  rigorous  proof  of  the  stability  of  a  mass-spring-dashpot  type  controller. 

2  Dissipative  System  Theory 

Intuitively,  a  dissipative  system  is  one  which  dissipates  energy.  In  other  words,  the  energy 
stored  in  such  a  system  will  be  less  than  the  energy  supplied  to  it  This  concept  can  be 
written  formally  in  terms  of  a  storage  function  which  generalizes  the  concept  of  energy.  The 
theory  in  this  section  is  taken  from  [5,  8,  9]. 
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Definition  2.1:  A  dynamical  system  E  is  defined  as  follows: 


(i)  U  is  the  input  space  and  consists  of  a  class  of  [/-valued  functions  on  the  positive  real 
line. 

(ii)  y  is  the  output  space  and  consists  of  a  class  of  K-valued  functions  on  the  positive  real 
line. 

(iii)  Both  Y  and  U  are  finite-dimensional  inner  product  spaces  and  both  y  and  U  are  closed 
under  the  shift  operator  i.e.  if  u{t)  6  U  then  u{t  +  T)  &U. 

(iv)  Define  R 2  :=  {(f2,ti)  G  RxR\  t2  >  tj}.  The  state  transition  function  <j>  :  R2  xX  xli 

X  and  defines  the  state  through  the  relation  x(t2)  —  (f>(t2,t1,x{ti),u)  This  function 
satisfies  the  usual  axioms  for  autonomous  dynamical  systems: 


(a)  Consistency  of  initial  condition:  f(t,t,x0,u )  =  x0 

(b)  Semigroup  property:  (f>(t2,  <1}  <f>(ti,  t0,  x0,  u),  u)  —  <f>(t2,  t0,  x0,  u) 

(c)  Causality:  u\ (t)  =  u2(t)  for  t0  <  t  <  t\  implies  (f}(t,to,xo,Ui)  =  (f)(t,tQ,xo,u2)  for 
t0  <t  <  fj. 

(d)  Time  Invariance:  f{tx  +  T, t2  +  T,x0,u *)  =  <f}(ti,t2,x0,u2)  for  all  T  >  0,t2  > 

<i,«a(0  =  +  T). 

(v)  The  function  r  :  X  x  U  >-»  Y  is  the  output  function:  y{t)  =  r(x(t),u(t)). 

In  view  of  the  system  time-invariance  (4d)  we  will  henceforth  use  to  =  0. 

Definition  2.2:  A  dynamical  system  is  said  to  be  dissipative  if  there  exists  a 
nonnegative  function  S  :  X  »->  R+,  called  the  storage  Junction  and  a  supply  rate  w  :  U  xY  •-» 
R  such  that  ,  for  all  t\  >  0, 

S(x0)  +  /  w(u,y)dt>S(x i)  (2) 

Jo 

where  u  G  U,xl  —  <f>(ti,  0,xo,u)  and  y  —  r(x,u). 

Note  that  the  storage  function  is,  in  general,  not  unique.  The  following  theorem 
shows  that  if  an  appropriate  supply  function  is  found,  it  is  not  necessary  to  actually  define 
a  storage  function. 


Theorem  2.1  [8]  If  there  exists  a  supply  function  w(u,y)  for  E  su  n  that 

rh 

/  w(u,y)dt  >  0 
Jo 

for  all  t\  >  0  then  there  exists  a  storage  function  S  such  that  (2)  is  satisfied  and  so  E  is  a 
dissipative  system. 


Moylan  and  Hill  [5]  define  a  useful  class  of  dissipative  systems  . 

Definition  2.3:  [5]  Given  P,  Q,  R  of  appropriate  dimension  with  Q,  R  symmetric, 
define  the  supply  rate 

Hu>  y)  =  (y.  pv)  +  ( y ,  Qu)  +  (u,  Ru)  (3) 

A  system  is  ( P ,  Q,  P)-dissipative  if 

ft  i 

/  w(u,y)  >  0  Vti  >  0. 

Jo 

Passivity,  or  positivity,  can  now  be  seen  to  be  the  special  case  where  P  =  R  =  0 
and  Q  —  I.  Strict  passivity  is  obtained  by  defining  R  =  —62I  and  Q  =  I.  Similiarly,  we 
obtain  input-output  stable  systems  by  choosing  P  =  —I  and  R  —  k2I  where  k  is  the  gain. 

Suppose  G  is  a  given  system,  for  which  we  wish  to  design  a  controller  H ,  arranged 
as  shown  in  Figure  1. 

The  feedback  system,  or  alternatively  the  pair  ( G,H ),  is  said  to  be  externally 
stable  if  u\,u2  G  £2(0,00,  U)  implys  y\,yi  €  £2(0,  00,  K),  and  there  is  a  maximum  ratio,  the 
£2  gain,  between  the  norm  of  the  input  and  the  norm  of  the  output. 

The  following  theorem  provides  a  simple  test  for  external  stability  of  intercon¬ 
nected  (P,Q,R)- dissipative  systems  . 


Theorem  2.2  [5][Theorem  1]  Consider  systems  G  and  H  connected  in  the  familar  feedback 
configuration  shown  in  Figure  1,  and  assume  that  H  is  (P2,Q2,  Ri)-dissipative  and  that  G 
is  (Px,  Qi,  Ri)-dissipative.  The  closed  loop  system  is  externally  stable  if 


a  _  Pi  +  P2  Q'x  ~  Q2 

v  ■“  Q'2  -  Qi  Ri  +  Pi  . 


(4) 


is  negative-definite.  □ 


The  well-known  passivity  and  small-gain  theorems  [3]  can  easily  be  derived  as  special  cases 
of  Theorem  2.2. 

In  order  to  study  stability  of  feedback  systems  with  more  general  supply  rates, 
we  state  an  observability  assumption. 

Assumption  1:  There  exists  some  T  >  0  and  a  non-negative  continuous 

function  a  :  R  >-*  R,  with  a(0 )  =  0  and  a(a)  >  0  for  a  >  0,  such  that  for  identically  zero 
input  and  any  initial  state  xo  ,  we  have 

/  y\t)y(t)dt>a( ||  i0||).  □ 

(For  finite-dimensional  linear  time-invariant  systems,  this  is  equivalent  to  the  standard  def¬ 
inition  of  observability.) 
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Theorem  2.3  Suppose  that 


(a)  A  given  system  G  has  a  supply  rate  wx  :  =  ( y ,  Pxy)  +  wc(u,y )  +  ( u,R\u )  with 

rh 

/  w1(u,y)dt>0, 

JO 


(b)  H  is  a  system  with  supply  rate  w2  :=  (y,  P2y )  +  wc(y,u )  +  ( u ,  R2u)  and 

rh 

/  w2(u,y)dt  >  0. 

Jo 

If  Ri  P\  and  P2  +  Rx  are  negative  definite  ,  and  both  G  and  H  satisfy  Assumption  1,  then 
the  origin  is  an  asymptotically  stable  equilibrium  point  of  the  closed  loop  system  (Figure  1 ) 
with  zero  external  inputs. 

Proof:  We  will  demonstrate  a  Lyapunov  function  for  the  closed  loop  system.  Let 
S i  and  S2  be  storage  functions  for  G  and  H  respectively.  Define  r  =  (r^rj)  and  similarly  y. 
We  have  ux  :=  rx  —  y2  and  u2  :=  r2  +  yx.  It  is  clear  that  the  interconnected  system  (Figure 
1)  is  dissipative  with  storage  function  S  :=  Sx  -f  S2  and  supply  rate 


w(r,y)  =  wx(u  1,1/! )  +w2(u2,y2). 

We  will  show  that  S  is  in  fact  a  Lyapunov  function  for  the  closed  loop  system. 

Since 

S{xjo)A  [  w(r,y)dt  >  S{xx)  (5) 

J  0 

it  follows  that  for  rx  —  r2  =  0, 

S  <  w  i(ui,  y  1)  +  w2(u2,y2) 

=  ^i(— 3/2, 3/1)  +  w2(yi,  y2) 

—  y'2^  1  +  -^2)2/2  +  y[  {P 1  +  R?)y  1 

<  0 

since  the  terms  arising  from  the  cross  rate  wc  cancel.  Thus  ,  S  is  negative  definite  (Assump¬ 
tion  1). 

It  now  remains  only  to  show  that  S  is  a  positive  definite  function. 

S(xq)  >  S(xi)  -  [  w(r,y)dt>  [  -y2 ( Rx  +  P2)y2dt  +  /  -y[{P\  +  Ri)yidt. 

Jo  —  Jo  Jo 

Since  both  plant  and  controller  satisfy  Assumption  1,  there  exists  a  positive  definite  function 
13  with 

S(xq)  >  0(||  Soil) 

and  so  S  is  a  Lyapunov  function  for  the  closed  loop  system.  The  result  follows.  □ 

Theorems  2.2  and  2.3  will  be  used  to  derive  controllers  for  uncertain  structures. 
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3  Rate  Sensors 


Suppose  we  have  a  structural  control  system  (1)  with  only  rate  sensors  i.e.  Cd  =  Ca  =  0. 
It  is  known  that  if  C'v  —  F  then  the  system  transfer  function  is  positive,  or  equivalently, 
the  system  is  (0, /,  0)-dissipative  [1].  The  passivity  theorem  can  be  used.  In  particular,  any 
strictly  positive  stable  controller  will  lead  to  a  stable  closed  loop  system.  It  is  not  necessary 
to  determine  the  model  (1)  beyond  ensuring  that  C'  =  F  and  that  the  dynamics  are  second 
order. 

In  this  section  we  extend  this  result  to  more  general  configurations. 

Theorem  3.1  Consider  (1)  with  Cd  —  Ca  =  0  and  suppose  that  there  exists  an  operator 
Q  :  U  ►— >  Y  such  that 

c;q  =  f.  (6) 

Then,  for  some  p2  >  0,  the  system  (1)  is  (— p2 /,  (5,0)  —  dissipative. 

Proof:  For  arbitrary  p  >  0,  define  the  supply  rate 
w(u,y)  :=  {y,-p2y)  +  (y,Qu). 


Define 


uT(t)  :=  PMt )  :=  {  “(i)-  ‘  ^  l 


and  let  yr{t)  be  the  output  which  corresponds  to  the  input  ur(t). 

r T  rT 

/  w(u,y)dt  =  /  -p2{yT,yr)  +  {yT,QuT)dt 

Jo  Jo 

r  oo  roo 

=  /  -p2{yT,yT)  +  (yT,QuT)dt+  p2(yT,yT)dt 

Jo  JT 

roo 

>  /  -p2(yT,yr)  +  (yT,QuT)dt 

Jo 

1  c°° 

=  7T  -p2{yT,iiT)  +  {yT,Qur)du 

Z/K  J — oo 

where  y  denotes  the  Laplace  transform  of  y. 

Defining  v(ju)  :=  ]<w(K  —  Mw2  +  D juj)-1  Fuj( ju>),  so  y  =  Cvv, 

[  wG(u,y)dt  >  J-  f  -p2(Cvv,Cvv)  +  (v,  Fur)duj 

J 0  Z7C  J  —  oo 

=  7T  r  v(ju)*{-p2CJCv  +  ~(K  -  Muj2  + 

Z7T  J —  oo  JUJ 

1  r°° 

=  2tt  /  v0uY{-P2c<i'Cd  + D)}v(]u>)duj. 
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If  we  choose  p2  small  enough  so  that  p2Cv'Cv  <  D ,  the  expression  in  square 
brackets  is  positive  definite  for  all  lj.  □ 

The  following  stability  theorem  is  now  immediate. 

Theorem  3.2  Assume  we  have  a  second  order  system  with  rate  sensors,  and  C'VQ  —  F,  as 
in  the  statement  of  the  previous  theorem.  If  a  given  controller  H  is  (1)  externally  Li-stable 
and  (2)  QH  is  a  positive  transfer  function,  then  the  closed  loop  system  is  stable. 

Proof:  The  system  is  (—  p2/,  Q,  0)-dissipative  for  some  p  >  0  (Theorem  3.1).  Let  k  be  the 
controller  gain  and  choose  a,  b  so  k  <  b2/a2  and  b 2  <  p2.  We  have  that  the  controller  is 
(—a2 1 ,Q' ,b2 1)- dissipative.  The  result  now  follows  from  Theorem  2.2.  □ 

This  result  is  more  general  in  several  ways  than  previous  results  depending  on 
positivity  of  a  structure.  The  plant  is  not  constrained  to  be  square  and  the  operator  Q  is 
not  required  to  bear  any  relation  to  a  Lyapunov  function  for  the  system.  However,  the  rate 
sensors  and  force  actuators  do  need  to  be  sufficiently  “collocated”  so  that  C'VQ  =  F  for  some 

Q- 

4  Structures  with  Displacement  Sensors 

In  many  applications,  sensing  is  performed  with  displacement  sensors.  If  the  control  is 
still  input  as  a  force,  then  the  system  transfer  function  is  by  no  means  positive,  even  the 
extended  sense  of  (0,£?,0)-  dissipativeness  discussed  in  the  previous  section.  We  would  like 
to  obtain  a  robustness  result  similar  to  that  obtained  in  the  previous  section.  One  approach 
to  this  problem  is  to  use  multipliers  (Figure  2).  Figure  2  can  be  shown  to  be  equivalent  to 
Figure  1.  The  multiplier  function  L  is  chosen  so  that  a  standard  result  such  as  the  passivity 
theorem  may  be  used  for  the  transformed  systems.  Unfortunately,  for  a  given  plant,  there 
is  no  guarantee  that  a  suitable  multiplier  function  exists,  nor  any  clue  of  how  to  choose 
the  multiplier  function.  Furthermore  difficulties  associated  with  ensuring  causality  of  the 
transformed  systems  may  arise.  Details  can  be  found  in  [3]. 

Looking  at  the  situation  from  a  physical  viewpoint,  replacing  rate  sensors  on  a 
structure  with  displacement  sensors  does  not  change  the  internal  dynamics  of  the  structure. 
Intuition  tells  us  that  the  system  should  remain  dissipative.  If  suitable  supply  and  storage 
functions  do  exist,  so  that  the  system  (1  )  with  displacement  sensors  is  dissipative,  Theorem 
2.3  can  then  be  applied  to  design  robust  controllers. 

In  order  to  keep  the  discussion  simple,  we  consider  the  situation  where  only 
measurements  of  displacement  are  made  i.e.  Cv  =  Ca  =  0.  We  note  first  that  the  definition 
of  supply  functions  can  be  extended  to  include  functions  which  involve  derivatives  of  the 
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inputs  and  outputs.  The  input  and  output  spaces  have  to  be  defined  appropriately  so  that 
the  supply  function  is  well-defined. 

Theorem  4.1  Consider  the  second-order  system  (1)  with  Ca  =  Cv  =  0,  assume  that  there 
exists  Q  so  C'dQ  =  F .  Let  c2  >  0  be  a  lower  bound  on  system  damping  so  that  c2M  <  D. 
Define  the  cross  rate 

wc(u,y)  :=  c2(y,Qu)  +  (y,Qu).  (8) 

Then  there  is  a  number  p2  >  0  so  that  this  system  is  dissipative  with  respect  to  the  supply 
rate 

wG(u,y)  =  - p2(y,y )  +  wc(u,y). 


Proof:  First  note  that  y  exists  for  all  outputs  y  so  that  wG  is  a  well-defined  supply 
function.  As  in  Theorem  3.1  and  in  [3]  we  will  show  that 

fT 

/  wG{u,y)dt  >  0 
Jo 

for  arbitrary  T  >  0  by  showing  that  its  transfer  function  is  positive.  Dissipativeness  of  the 
system  with  this  supply  rate  will  then  follow  from  Theorem  2.1. 

Defining  ur{t)  :=  Ptu(1)  as  in  equation  (7),  let  yr(t)  be  the  output  which  corre¬ 
sponds  to  the  input  uT(t). 

rT  fT  2 

/  wG(u,y)dt  =  /  -p(yT,yT)  +  c(yT,QuT)  +  (yT,QuT)dt 

Jo  Jo 

r°°  f°° 

—  /  —  P2{yT,VT)  +  c2{yT,QaT)  -I-  {yT,QuT)dt  -f  /  p  (yT,yT)dt 

Jo  Jt 

>  /  - p2{yr,yT )  +  C2(yr,  Qut )  -I-  (yT,  QuT)dt 

Jo 

1  f°° 

=  7T  -p\yf,yf)  +  c2(yT,QuT)  +  (syT,QuT)duj. 

Z  7T  J — oo 

Defining  v(ju>)  :=  ( K  —  Moj2  +  Fut  'joj), 

f  wG(u,y)dt  >  ~  [  v(]U})*[-p2Cd'Cd  +  c2(K  -  Muj2)  +  Doj2)v{]oj)dui . 

Jo  In  J -oo 

If  we  choose  p2  so  that  p2Cd'Cd  <  c2K,  the  expression  in  square  brackets  is 
positive  definite  for  all  u>  and  for  arbitrary  T, 


rT 

/  wG(u,y)dt  >  0 
Jo 


as  required.  □ 


Theorem  4.2  Consider  the  structural  system  defined  in  the  previous  theorem.  Assume  ob¬ 
servability.  The  closed  loop  system  (Figure  1)  will  be  stable  if  a  linear  controller  is  (a)  stable 
(b)  dissipative  with  respect  to  the  supply  rate 

w(u,  y)  =  wc(y,  u )  =  c2(u,  Qy )  +  ( u ,  Qy )  (9) 

where  Q,c2  and  wc  arc  as  defined  in  the  statement  of  the  previous  theorem.  The  input  space 
U  is  restricted  to  inputs  with  u  6  C2{U). 

Proof:  Define  the  positive  number  p2  as  in  the  previous  theorem.  Since  the  controller  is 
stable,  it  is  (— a2 1, 0,  62/)-dissipative  where  we  choose  b 2  <  p2  and  b2/a2  large  than  the 
controller  gain.  Choose  an  observable  realization  for  the  controller.  If  condition  (b)  also 
holds,  then  it  follows  that  the  origin  is  a  stable  equilibrium  point  of  the  closed  loop  system 
with  U\  —  U2  —  0(Theorem  2.3  ).  Since  both  plant  and  controller  are  linear  and  finite 
dimensional,  it  follows  that  the  closed  loop  system  is  also  externally  stable.  □ 

Errors  in  the  model  of  the  system  stiffness  K  or  actual  damping  larger  than  the 
estimate  will  not  destabilize  the  closed  loop  system.  Other  than  ensuring  C'dQ  =  F,  the  only 
knowledge  of  the  structure  required  in  designing  a  stabilizing  controller  is  a  lower  bound  on 
the  damping  so  that  the  controller  satisfies  (9).  If  the  actual  system  damping  is  greater  than 
the  estimate,  the  closed  loop  system  will  still  be  stable. 

Note  that  this  result  could  also  have  been  obtained  through  defining  the  multiplier 
function  L  —  Q(c2-\-s).  We  feel  however  that  there  are  several  advantages  to  using  dissipative 
system  theory.  The  first  one  is  that  supply  rates  such  as  (u,y)  may  be  defined  without 
the  problems  of  causality  and  stability  associated  with  the  multiplier  function  s.  More 
importantly,  we  think  that  the  the  choice  of  supply  rate  is  clearer  when  it  is  known  that  the 
system  “dissipates  energy”  in  some  sense. 

5  Example:  Second  Order  Controllers 

In  order  to  illustrate  the  ideas  discussed  above,  consider  the  second  order  dynamic  system 
with  displacement  sensors  and  C'dQ  —  F  for  some  Q  as  discussed  in  the  previous  section. 

Suppose  we  wish  to  control  this  system  with  a  combination  of  direct  feedback 
and  the  displacements  of  a  second  order  system.  Mathematically,  the  proposed  controller 
has  the  form 


M2x2{t )  +  D2x2(t )  +  K2x(t)  =  F2u2(t) 
y2(t)  =  Cj2x2(t)  +  Eu2(t). 


(10) 


Since  both  plant  and  controller  are  stable  some  readers  might  expect  stability  of 
the  closed  loop  system  to  follow  if  both  plant  and  controller  have  positive  damping.  This  is 
not  true,  as  the  following  simple  example  illustrates. 


Example  :  x{t)  +  8x(f)  +  100i(t)  =  uj(t)  x2(f)  -f-  .001x2(t)  +  .lx2(t)  =  u2(t) 


Ui  —  —  X2 


U2  =  V 2  ~h  X 


While  for  damping  d  —  8,  the  closed  loop  is  stable,  increasing  the  plant  damping  to  12  will 
lead  to  an  unstable  closed  loop. 

However,  by  using  the  results  from  the  previous  section  we  can  obtain  a  controller 
(10)  which  is  unconditionally  stable. 

Theorem  5.1  Suppose  we  have  a  observable  second  order  plant  (1)  with  Cv  —  Ca  =  0 
and  (2)  there  exists  Q  such  that  C'dQ  ~  F  .  If  the  controller  (10)  is  implemented  so  that 
Cd2Q'  —  ~E2  and  the  bias  E  satisfies  QE  >  F2'X'F2  where  X  >  K2~l ,  then  the  closed  loop 
system  is  stable. 

Proof:  From  Theorems  4.1  and  4.2  the  result  will  follow  if  we  can  show  that  the 
given  system  is  dissipative  with  respect  to  the  supply  rate  (9)  for  arbitrarily  small  c2  >  0. 

Define  v(ju>)  :=  (  —  M2uj2  +  D2j u  +  K2)~1ut{]oj)^  and  set  E'Q'  =  F2XF2  -f  E2 
where  X,  E2  are  as  yet  undetermined.  Using  the  same  frequency  domain  technique  as  in 
Theorems  3.1  and  4.1,  we  obtain 


[T  roo 

/  wc(u,y)dt  >c2  v*  (  jw)Z(ju))v(ju >)  +  u?*  {ju)G 2ut(jw)  du 

JO  J —oo 


where 


Z(juj)  =  [K2XK2  -  K2]  +  oj2{D2/c2  +  M2-  M2XK2  -  K2XM2  +  D2XD2\  +  w*M2XM2  . 

Choose  c2  small  enough  that  (1)  c2M  <  D  and  (2)  the  coefficient  of  u>2  is  positive  definite. 
If  XK2  >  /,  then  Z  is  positive,  and  the  controller  is  dissipative  with  supply  rate  wc  for  all 
E2  >  0.  Therefore,  the  closed  loop  is  stable  (Theorem  2.3).  That  is,  the  closed  loop  system 
is  stable  for  all  E  which  satisfy 

QE  >  F2K2~lF2  .□  (11) 

Does  this  design  have  any  physical  meaning?  The  answer  is  positive.  Consider 
the  special  case  where  the  controller  and  the  system  have  the  same  number  of  states,  and 
furthermore  assume  that  all  states  are  measurable  and  the  actuators  are  collocated  with  the 


sensors. 


Mxi(t)  +  Dxi(t)  +  Kxi(t )  =  ui  =  rx(t)  -  y2(fi 
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Adjust  the  controller  gain  so  F2  =  K2: 


M2x2(t)  +  D2x2(t)  +  K2x2(t)  —  K2u2  —  K2{r2(t)  +  yi(t)) 
y2  =  -K2x2  +  Eu2  =  K2x2  +  E(yi  +  r2) 


(13) 


Defining  x  :  =  (x1,i2)  we  can  write  the  closed  loop  system  as 

K  +  E  -if2  1 


M  0 
0  M2 


x(t)  + 


D  0 

o  d2 


x{t)  + 


k2 


k2 


x(t)  = 


1  ~E 

o  k2 


t{t). 


This  is  a  stable  second  order  system  as  long  as  the  matrix 

Kt  = 


K  +  E 
K2 


-K  2 
K2 


is  positive  definite  i.e.  it  represents  the  stiffness  matrix  for  two  “springs”  K  and  K2  connected 
in  series.  This  is  true  if  E  >  K2.  Noting  that  F2  —  K2  and  Q  —  / ,  this  inequality  is  identical 
to  the  condition  (11)  arrived  at  using  the  dissipativeness  criterion. 

If  we  have  scalar  systems,  the  choice  E  —  k2  leads  to  a  stable  closed  loop  system. 
Setting  the  reference  inputs  to  zero, 


tti  =  ~V2  =  k2(x2  -  xi) 


The  controller  in  this  case  reduces  to  a  spring-mass-dashpot  system  connected  in  series  with 
the  system  mass  as  shown  in  Figure  3. 

For  d2  =  0,  the  system  is  just  two  spring-masses  connected  in  series.  As  k2  is 
increased,  the  stiffness  of  the  system  increases.  For  d2  >  0,  the  system  is  always  stable  for 
k2  >  0,  and  increased  amounts  of  energy  are  dissipated  in  the  dashpot  as  d2  is  increased.  By 
adjusting  the  control  system  parameters,  k2,d2,m2  ,  we  can  design  the  closed  system  to  be 
under-,  critically  or  over-damped. 

6  Conclusions 

We  have  used  results  in  dissipative  system  theory  and  Lynpunov’s  Second  Method  to  develop 
stability  theorems  which  generalize  the  passivity  and  small-gain  theorems.  Simple  conditions 
for  robust,  model-independent  controllers  which  stabilize  non-passive  structures  were  derived. 
Throughout  the  paper  we  have  used  the  second-order  form  of  the  dynamical  equations.  Only 
displacement  and  velocity  feedback  were  studied.  However,  the  extension  to  accelerometers, 
and  to  combinations  of  these  three  types  of  sensors  is  straightforward. 


Although  for  displacement  feedback,  a  lower  bound  on  the  system  damping  is 
required,  the  controller  is  otherwise  independent  of  the  system  parameters.  It  should  be 
emphasized  that  this  is  a  robustness  result  with  respect  to  structural  uncertainty  in  the 
absence  of  measurement  uncertainty  and  other  contributing  factors. 

Control  performance  is,  of  course,  dependent  on  the  system  characteristics.  Knowl¬ 
edge  of  the  system  model  can  always  help  improve  a  controller  design.  Future  work  will  ex¬ 
plore  combining  satisfaction  of  the  appropriate  dissipative  condition  with  design  techniques 
such  as  LQG  regulator  and  H methods.  Finally,  the  controller  has  been  formulated  from 
the  continuous-time  setting.  Actual  implementation  of  the  controller,  however,  most  likely 
requires  usage  of  a  digital  computer.  In  future  work,  effects  of  sampling  and  time  delays 
will  be  addressed.  Other  practical  issue  that  can  also  affect  the  control  performance  such  as 
measurement  noises,  and  actuator  and  sensor  saturation  limits  will  be  investigated. 
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Figure  2:  Feedback  System  with  Multipliers 


SYSTEM  CONTROLLER 

( A  Real  Spring-Mass  System)  (A  Virtual  Spring-Mass-Dashpot  System) 

f - T - I 


I _ 1 _ I 

Figure  3:  A  simple  spring-mass  system  with  a  simple  spring-mass  controller 
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